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We have performed precision microwave spectroscopy on ultra-cold *^Rb confined in a magnetic 
trap, both above and below the Bose-condensation transition. The cold collision frequency shifts 
for both normal and condensed clouds were measured, which allowed the intra- and inter-state 
density correlations (characterized by sometimes controversial "factors of 2") to be determined. 
Additionally, temporal coherence of the normal cloud was studied, and the importance of mean- 
field and velocity-changing collisions in preserving coherence is discussed. 
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With the advent of modern cooling and trapping tech- 
niques, the cold collision regime has become readily ac- 
cessible. In the cold collision regime quantum statistical 
effects due to particle indistinguishability dominate scat- 
tering processes. The symmetrization requirement for 
identical particles in an ultra-cold Bose gas enhances the 
probability of finding two particles nearby, causing den- 
sity fluctuations. At lower temperatures the statistical 
nature of the Bose gas causes the phenomena of Bose- 
Einstein condensation, where all atoms in the conden- 
sate share the same wavef unction, suppressing the den- 
sity fluctuations found in a noncondensed sample. 

Suppression of second-order density fluctuations in a 
condensate has been measured through analysis of the 
expansion energy of condensates [0, |^] . In a separate ex- 
periment, the suppression of third-order density fluctua- 
tions was probed by comparing the three-body loss rate 
of a condensate to that of a normal cloud |^ . The effect 
of cold collisions has also been measured as a density- 
dependent energy shift in atomic fountain clocks jH ||, || . 
These shifts are quite small (~0.1-10 mHz) due to the low 
densities at which the clocks operate, but are measurable 
because of their high precision. The uncertainty associ- 
ated with these collisional shifts can be problematic; in 
fact the next generation of atomic fountain clocks are 
based on *^Rb rather than ^^'^Cs because the collisional 
shift of ^^Rb is ~ 30 times smaller. In recent ultra-cold 
hydrogen experiments the cold collision shift provided 
the signature of Bose-condensation; below the transition 
a large frequency shift of the 1S-2S transition was seen, 
reflecting the high density of the condensate Q . 

In this paper we report precision microwave spec- 
troscopy performed on ultra-cold and condensed ^^Rb 
atoms confined in a magnetic trap. Due to the high den- 
sities achievable in a magnetic trap, the collisional energy 
shifts were 10^ greater than those in ®^Rb atomic clocks, 
allowing a high-precision measurement of the shifts of the 
magnetically trappable states to be made with relative 
ease. The collisional shifts for both a normal and con- 
densed sample were measured, providing a useful probe 
of the quantum statistics of the system. Additionally, 
magnetic confinement permits long interrogation times. 



allowing us to characterize temporal coherence of the nor- 
mal cloud under various experimental conditions. Com- 
parison of measured coherence times with a coUisionless 
numerical simulation suggests that collisions preserve co- 
herence in normal clouds. 

Precision spectroscopy of trapped samples is difficult 
because atom trapping relies on spatial inhomogencity 
of the atomic energy levels. Spatial inhomogencity of the 
energy levels broadens the transition frequency, thus lim- 
iting the precision attainable through spectroscopy. In 
other work, this difficulty has been avoided by confining 
atoms in a blue-detuned optical dipole trap, where atoms 
spend little time interacting with the trapping fields 
In this work spatial inhomogencity of the transition fre- 
quency was minimized through the use of a pair of energy 
levels which experience the same trapping potential. At 
a magnetic field of ^ 3.23 G, the |1) = |F = 1, to/ = —1) 
and |2) = |_F 2,r7i/ ~ 1) hyperfine levels of the 5Si/2 
ground state of ®^Rb experience the same first-order Zee- 
man shift. For a normal cloud at 500 nK, each energy 
level is Zeeman shifted by ~10 kHz across the extent of 
the cloud; however, at 3.23 G the differential shift of the 
two levels across the cloud is ~1 Hz. Compared to the 
differential Zeeman shift, the energy shift due to cold col- 
lisions is then a relatively large effect at high densities, 
making measurements of collisional shifts in this system 
possible. The small inhomogeneity allows for long coher- 
ence times, ~2 seconds and longer for low-density clouds, 
making this system attractive for precision measurements 
as well as for the study of condensate coherence in the 
presence of a thermal cloud. 

The experimental setup has been previously described 
1^ and will be briefly summarized here. Approxi- 
mately 10^ ^^Rb atoms are loaded into a vapor cell 
magneto-optical trap (MOT). The atoms are then op- 
tically pumped into the \F — 1) state by turning off 
the repump beam while MOT beams remain on. Then 
the trapping beams are turned off and the MOT coils 
are ramped to a high current forming a 250 G/cm gra- 
dient to trap 1 1,-1) atoms in the quadrupole field of 
the coils. The quadrupole coils are mounted on a lin- 
ear servo-motor controlled track which then moves the 
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coils 44 cm, from the MOT region to a lofFe-Pritchard 
trap in the ultra-high vacuum region of the system. The 
loffe-Pritchard trap consists of two permanent magnets 
which provide a 450 G/cm radial gradient. Two pairs of 
electromagnetic coils, a pinch and a bias, provide con- 
finement in the axial direction, which is aligned perpen- 
dicular with respect to gravity. At a typical bias field of 
3.23 G, atoms in the |1, —1) state experience {230, 230, 7} 
Hz trap frequencies. The sample is further cooled by rf 
evaporation, and condensates of up to 10^ atoms can be 
formed. Imaging is performed by the use of adiabatic 
rapid passage to transfer atoms from the |1, — 1) state to 
the |2,— 2) state. Anti-trapped |2,— 2) atoms rapidly ex- 
pand for 2-5 ms and then are imaged through absorption 
by a 20 ^s pulse of 5Si/2 |2, -2) ^ 5P3/2 |3, -3) light. 

A two-photon microwave-rf transition is used to trans- 
fer atoms between the |1) and |2) states. A detuning 
of 0.7 MHz from the |2,0) intermediate state provides 
a two-photon Rabi frequency of 2.5 kHz. Ramsey 
spectroscopy of the |1) — s- |2) transition is performed by 
measuring the total number of atoms remaining in state 
|1) after a pair of ^ pulses separated by a variable time 
delay are applied The frequency of the resulting 

Ramsey fringes is the difference between the transition 
frequency 1^12 and the two-photon drive frequency. In 
previous work we measured local variations of 1^12 by de- 
tecting the number of atoms remaining in state |1) at 
specific spatial locations along the axis of the normal 
cloud |9|. By analyzing the spatio-temporal variations 
of 1^12, combined with the measured evolution of the |1) 
state after a single ^ pulse, we were able to spatially 
resolve the evolution of spin waves [ pT| . In this work, 
in order to perform measurements of 1/12 insensitive to 
spin waves, one of the following two techniques was used. 
With one technique the entire cloud, rather than specific 
spatial locations, was monitored to average out the ef- 
fects of spin waves. Alternatively, Ramsey spectroscopy 
was restricted to interrogation times short compared to 
the spin wave frequency [ p^ . 

One effect which shifts the transition frequency 1^12 is 
the differential Zeeman shift. The Breit-Rabi formula 
predicts a minimum in 1/12 at Bq — 3.228917(3) Gauss, 
thus the differential Zeeman shift between the |1) and |2) 
energy levels is first-order magnetic field independent at 
B = Bq. The differential Zeeman shift about Bq can be 
approximated as 1^12 = Vmin+PiB—B^Y Measuring 
1^12 for different magnetic fields allows us to calibrate our 
magnetic field from the expected dependence, see Fig. ||. 
By working in the vicinity oi B = Bq we greatly reduce 
spatial inhomogeneity of 1/12 and also become first-order 
insensitive to temporal magnetic field fluctuations. 

A second effect which shifts z/12 arises from atom-atom 
interactions. In the s-wave regime, where the thermal 
de Broglie wavelength of the atoms is greater than their 
scattering length, atoms experience an energy shift equal 
to a^^^an, where a is the two-particle correlation at 
zero separation, n is atom number density, a is the scat- 
tering length, and m is the atom mass. Therefore for a 
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FIG. 1: Differential Zeeman shift at low magnetic fields for the 
jl) ™» |2) {\F = l,m/ = -1) ^ \F = 2,mf = 1)) transition. 
The solid line is the predicted splitting from the Breit-Rabi 
formula. The inset plot expands the bias field region where 
most studies are performed. 



two-component sample the expected energy shift of each 
state is 



d/ii = (aiiaiini + 012012^2) (1) 
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where ni and n2 are the |1) and |2) state density and 
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The shift of the transition frequency in Hz can then be 
written as 



Az/12 = 



2h, 
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m 

h . 

= — n(a22a22 - aiiaii + 



(2ai2ai2 - aiioii - 0:22022)/), 



(4) 



where / — "^""'^ and n — ni + n2. 

For noncondensed, indistinguishable bosons, a ~ 2 due 
to exchange symmetry, therefore a"f = aJI = 2 in a cold 
normal cloud (where the superscript c or nc refers to 
condensed of noncondensed atoms respectively). Distin- 
guishable particles do not maintain exchange symmetry, 
making q;"| = 1 for an incoherent two-component mix- 
ture. However if a two-component sample is prepared by 
coherently transferring atoms from a single component. 
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such as in Ramsey spectroscopy, then the excitation pro- 
cess maintains exchange symmetry, and we might expect 
a f = 2 In this scenario the coUisional shift should 
be calculated using a"{^ = ~ '^12 — 2, leading to a 
predicted frequency shift of 



Aj/12 = — n(a22 - an + (2ai2 - an - a22)/). (5) 

TO 

This result can also be obtained by solving the trans- 
port equation From spectroscopic studies Q 
the three ^^Rb scattering lengths of interest have been 
determined to be a22 — 95.47ao, 012 — 98.09ao, and 
an = 100.44ao, where uq is the Bohr radius. The fre- 
quency shift can then be written as 



2/1 

Ai/12 = —aon(-4.97 4-0.27/). 

TO 



(6) 



If on the other hand the |1) and |2) states do not main- 
tain exchange symmetry, such that a"| — 1, then the 
frequency shift would instead be 



Ai/12 = —aon(-^4.97- 97.82/). 

TO 



(7) 



These two models are clearly distinguished by the depen- 
dence of V12 on /. 

When we perform Ramsey spectroscopy with a pair of 
^ pulses, the populations of the |1) and |2) states are 
equal, and thus f — during the interrogation time. 
From Eq. (4) it is apparent that with f — the coUi- 
sional shift is sensitive only to a"^ and an terms. For 
these measurements the bias field was set to Bq, and the 
transition frequency was measured for a range of den- 
sities. To adjust density of the sample, the number of 
atoms in the initial MOT load was varied. All normal 
cloud data was taken at the same temperature of 480 nK, 
and all condensate data was taken with high condensate 
fractions in order to minimize effects due to the normal 
cloud. The density for the normal cloud was found by fit- 
ting Gaussian profiles to absorption images of the clouds 
and extracting the number, temperature, and density. To 
measure condensate density, Thomas-Fermi profiles were 
fit to absorption images of the condensates and the total 
number, A'o, in the condensates and the Thomas- Fermi 
radius along the long axis, Z, were extracted. 

The results of this measurement are shown in Fig. |^. 
Comparing the coUisional shift measured for the normal 
cloud to that measured for a condensate gives ct^f/a'i^ = 
2.1(2). If instead we assume a^f = 2 and a^j = 1, then 
the data for both the condensate and normal cloud can 
be used to obtain a value for the difference in scattering 
lengths of a22 — an = — 4.92(28)ao, in agreement with 
values determined from molecular spectroscopy ]l7[ |. 

Many systematics can plague density measurements 
made through absorption imaging. In order to test in- 
dependently our density calibration for both the nor- 
mal and condensed samples, each of which can suffer 
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FIG. 2: Measurement of the cold collision shift. Solid and 
open circles represent measurements of the normal cloud and 
condensate respectively. The solid line is a fit to the normal 
cloud data = 0.1(0.4)-3.9(0.3)10"^^n; the dashed line is 
a fit to the condensate data Ai/12 = -0.1(1.4)-1.9(0.2)10"^^n 
where Aui2 is in Hz and n is in cm~"^. 



from different errors, we used the Bose-condensation phe- 
nomenon. The density of normal clouds was tested 
through measurement of the critical temperature, and 
condensate density was tested with the Thomas-Fermi 
approximation. Assuming that disagreements are due 
only to errors in estimation of atom number, the worst 
case scenario, leads us to reduce normal cloud density by 
11(4)% and increase condensate density by 11(3)%. Ad- 
justing the cold collision shifts accordingly would yield a 
worst-case corrected value of a^f /a'^^ = 1-7(2) [Q. The 
adjusted normal and condensate density shifts can be 
combined as above to give a value for the difference in 
scattering lengths of a22 — an — — 4.85(31)ao; not signif- 
icantly different from our unadjusted measurement. 

The remaining significant systematic is atom loss due 
to 1 2)- 1 2) coUisional dipolar relaxation. In order to min- 
imize effects of this loss, interrogation times were kept 
as short as possible. Nevertheless for the highest density 
condensate measurements the |2)-|2) loss causes the total 
density to drop by 3% in 20 ms, the maximum interroga- 
tion time. For all other densities the loss was no larger 
than this, and in most cases much smaller. Finally, the 
|1) and 1 2) states begin undergoing spatial separation in 
the condensate after the first | pulse however, the 
timescale for the separation is much longer than our 20 
ms interrogation time. 

Exchange symmetry between the |1) and |2) states can 
be tested by working at a fixed density and varying the 
relative |1) to |2) population by varying the length of the 
first Ramsey pulse |2^. In this case the first term in 
Eq. (4) will be constant and the measurement will test 
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FIG. 3: Testing the exchange symmetry between the |1) and 
|2) states. The transition frequency is measured as f — 
is varied for a normal cloud at fixed peak density of 7 x 10^^ 
cm~^ and temperature of 510 nK. The solid line is the fit, 
which yields ai2/Q!ii,22 = 1-01(2), which is to say, inter- and 
intra-state density correlations are quite accurately the same. 
The dotted line indicates the expected slope for a"! /^ii 22 = 
1/2. 



a"| and 012 as well as the a^f and an terms (see Eq. (6) 
and (7)). To minimize systematics the interrogation 
times were kept short, making precise frequency determi- 
nation difhcult. Nevertheless, our measurement (Fig. ||) 
indicates ai2/aii,22 = 1-01(2), where we have used the 
spectroscopically determined scattering lengths. This 
clearly indicates that exchange symmetry is maintained 
between the |1) and |2) states. A similar measurement 
was made on the \F — l,mf — 0) \F = 2, m/ = 0) 
transition by Fertig and Gibble 

As a thought experiment, imagine distinct thermal 
populations of |1) and |2) atoms, separately prepared, 
then mixed together, with the energy of interaction (pro- 
portional to a"|) measured for instance calorimetrically. 
Surely in this case the density fluctuations in state |1) 
and in state |2) would be uncorrelated, and a"| would 
be determined to be 1, not 2. We lack the experimental 
sensitivity to make such a calorimetric measurement, and 
our Ramsey-fringe method of measuring energy differ- 
ences obviously would not work for incoherent mixtures. 
We speculated, however, that if a"| = 2 for coherent 
superpositions, and if a"| = 1 for incoherent mixtures, 
then for partially decohered samples, a"| would take on 
some intermediate value. So by performing a measure- 
ment similar to that in Fig. |^ we might expect to see a 
more negative slope for a partially decohered sample; al- 
ternatively a frequency chirp in the Ramsey fringes may 
be seen as the sample decoheres. 

We probed the time evolution of a"| in a way similar 
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FIG. 4: A data set of Ramsey fringes probing for frequency 
shifts as a function of coherence. For this measurement nor- 
mal clouds at a temperature of 480 nK and a peak density of 
3.2 X 10^^ cm~^ were used. 

to Fig. ^; however rather than varying / we set / ~ 0.8 
then measured i>i2 with long interrogation times, looking 
for a frequency chirp as the fringe contrast decreased. 
This method has the advantage that there is a relatively 
small 1 2) state population, so effects arising from |2) loss 
are minimized. Seven data sets were taken for this mea- 
surement; an example is shown in Fig. ^. By allowing a 
linear frequency chirp in the fit of the Ramsey fringes, the 
frequency shift can be constrained to — 0.2(^ Hz by the 
time the fringe contrast has reduced to 1/e [£1| . However 
if we hypothesize that q;"| goes from 2 to 1 linearly as 
fringe contrast goes from 100% to 0% we would expect a 
frequency shift of —20(2) Hz as the fringe decayed, while 
the experimental limit is 40 times smaller. Clearly this 
appealing but unrigorous model is far too naive. 

Ramsey spectroscopy not only allows us to probe the 
energy difference between the two states, but also per- 
mits the measurement of the coherence between the two 
states. Coherence measurements were performed using 
the same time domain method used to measure V12 ; how- 
ever, interrogation times were extended until the Ramsey 
fringe contrast was lost. The resulting data was fit to a 
e-(*/T) decay, where r is the coherence time. Addition- 
ally, fitting allowed for both a loss in total atom number 
and a linear frequency chirp of 1^12. To ensure that the 
1/e atom loss times were much longer than the coherence 
times fractional transfers (/ ~ 0.8) to the |2) state were 
used for the high density data points. The results of co- 
herence measurements for different magnetic fields and 
at three different densities are shown in Fig. 

We expect that the primary source of decoherence is 
inhomogeneity in P12 across the cloud, to which both 
the differential Zeeman and coUisional shifts contribute. 
The collisional shift scales directly with density, and thus 
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FIG. 5: Coherence times in a normal cloud. The temperature 
was ^ 480 nK and the normal cloud peak density for each 



: (b) 










1 




\ 


li 


} 

i 


i 




i 










1 


2.9 


3.0 


3.1 


3.2 


3.3 


" 1 

3.4 







1 






i " 








2.9 


3.0 3.1 3.2 3.3 
Bias Field (Gauss) 


3.4 



plot is: (a) 4 x 10 



12 



cm 



(b) 1.3 X 10^^ 



cm 



(c) 3.2 X 10 
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cm The solid line in each plot corresponds to the 1/e times 
obtained from a numerical simulation (see text). The noise 
in the simulation data is due to the finite number of particles 
used and the random initialization. 



provides a Gaussian-shaped vi2 profile across the cloud 
Inhomogeneity arising from the differential Zeeman 
effect depends on the bias field. The magnitude of the 
magnetic field near the bottom of the trap can be writ- 
ten as B{z) = ~\- Bbias, thus the inhomogeneity due 

to the Zeeman shift is 1^12 (x -^4— -I- {Buas — Bq)B" z^. 
By setting Bi,ias{>, <, =)Bo the curvature of 1^12 can be 
adjusted to be positive, negative, or nearly zero, respec- 
tively. This allows inhomogeneity due to the coUisional 



shift to be roughly cancelled by an opposing Zeeman in- 
homogeneity [pl . This cancellation can be seen in Fig. 
as the density is increased from (a) to (c), the inhomo- 
geneity induced by the coUisional shift increases, so that 
a larger opposing Zeeman inhomogeneity is necessary for 
cancellation. Therefore the bias field for peak coherence 
time decreases as cloud density increases. 

In an attempt to compare the measured decoherence 
times to the known spatial inhomogeneity of the tran- 
sition frequency, we performed the following numerical 
simulation: consistent with a Maxwell-Boltzman distri- 
bution, we randomly assign initial positions and veloci- 
ties to ten thousand simulated atoms. Ignoring the effects 
of collisions, we calculate the three-dimensional trajec- 
tory of each atom for several simulated seconds, keeping 
track of the time integral of ^12 along the trajectory. At 
each point in time, we calculate the spatially integrated 
transverse magnetization and, as inhomogeneities cause 
this magnetization to wash out, find the time it takes 
the integrated transverse magnetization to reduce to 1/e 
of its original value This model should correctly 

account for the effect of motional averaging except that 
all coUisional effects are explicitly excluded. The result- 
ing modelled damping times are plotted as a solid line 
along with the experimentally measured damping times 
in Fig. Ij. While the model does a reasonable job predict- 
ing the value of the bias field for which the coherence time 
peaks, it consistently underestimates (in some cases by 
a factor of eight) the actual value of the coherence time. 
Our model neglects both velocity-changing collisions and 
the exchange- type collisions that lead to spin waves; it 
appears that these effects contribute significantly to pre- 
serving coherence across the trapped atom cloud. 

The extreme aspect ratio of our trapped cloud com- 
plicates a proper quantitative analysis of the effects of 
collisions on coherence. Along the axial direction, our 
previous work has shown that the effects of spin-waves 
are to keep local magnetization across the cloud from 
straying too far from its spatially averaged value [see in 
particular Fig. 2(a) of Ref. ^^]- In the radial directions, 
the motional oscillation frequency exceeds the mean-field 
exchange frequency in the cloud, and thus the effects 
of spin- waves are probably not relevant. On the other 
hand, velocity-changing collisions are likely important - 
in their absence, atoms with small transverse energies 
would stay near the axis of the trap, while atoms with 
large transverse energies would preferentially sample the 
larger magnetic fields and lower densities further from the 
axis. Velocity-changing collisions which re-randomize the 
transverse trajectories of these different classes of atoms 
(before they have a time to accumulate a radian or more 
of relative phase-difference) will serve to extend the co- 
herence time of the sample. In the highest density data 
set presented [Fig. ^(c)] the mean elastic collision rate 
was 74 Hz, which should be compared the measured co- 
herence times of around 0.5 to 1.2 seconds. 

It is interesting to consider the usefulness of magnet- 
ically trapped atoms for precision metrology. Peak co- 
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herence times of approximately 2.5 s were realized with 
cold, low density samples. An interrogation time of 2.5 
s provides a 0.2 Hz linewidth, which naturally leads one 
to consider using such a system for precision measure- 
ments. By working at B^ias = Bq, coherence times are 
slightly reduced. However, perturbations of 1^12 due to 
magnetic field fluctuations become very small, on the 
order of 4 mHz for current typical experimental condi- 
tions. With careful design of the confinement coils, the 
current supply, and magnetic shielding of external fields, 
it should be possible to suppress fluctuations of the bias 
field below 1 mG, reducing Zeeman-induced frequency 
shifts to below 0.1 mHz. Perturbations of 1^12 due to the 
cold collision shift are more significant; shot-to-shot den- 
sity fluctuations will introduce frequency noise, so it is 
advantageous to work with the minimum possible den- 
sity. However as atom number is reduced the maximum 
signal-to- noise ratio will decrease due to shot-noise , 
therefore the optimum strategy is to work with an atom 
number such that the frequency uncertainty due to shot- 
noise is on the order of the uncertainty due to density 
fluctuations. 

For example a normal cloud of 400 nK and a peak den- 
sity 1.5 X 10^^ cm~^ has 6 x 10'* atoms. The shot-noise- 
limited signal-to-noise ratio is then 245:1. With an inter- 
rogation time of 1 second the single-shot statistical uncer- 
tainty is then 0.65 mHz; including the effects of decoher- 
ence and atom loss will increase this to approximately 
0.9 mHz. Assuming that shot-to-shot number fluctua- 
tions are also shot-noise limited, then statistical uncer- 
tainty from the density shift is 0.84 mHz. Combining 
these gives a total single-shot uncertainty of 1.24 mHz; 
with our current 30 second evaporation time, the duty- 
cycle is such that an absolute precision of 6.8 mHz/\/Hz 
can be realized. This corresponds to a relative precision 
of 1 X 10^^^ /\/Hz, which in terms of measurement pre- 



cision does not reach the level of atomic fountain clocks. 
This system however has the advantage that small energy 
shifts can be measured in a compact, stationary spatial 
position; a 400 nK cloud occupies approximately only a 
1040 X 32 X 32 /im region of space. It is certainly feasible 
to perform spectroscopy at this level within 100 fim of 
a surface, which might allow the measurement of short 
range atom-surface interactions. 

In future work we plan on measuring coherence times 
in finite temperature condensates to study the role of the 
normal cloud in decoherence. We anticipate this system 
will be quite rich due to the existence of spin waves in the 
normal component, phase separation in the condensate, 
and the interaction between the two. The complexity 
of this interaction may partially account for the anoma- 
lous density shift of the condensate seen in the hydrogen 
experiments ^ p^ . 

We have demonstrated precise spectroscopy in an 
ultra-cold magnetically trapped gas. This permitted 
measurement of the cold collision shift in both a conden- 
sate and a normal cloud, allowing a probe of the quantum 
statistics of the system. Working at low densities mini- 
mizes the effect of the collisional shift, allowing long co- 
herence times and precise determination of (^12; however 
the measurement of any quantity not related to atom- 
atom interactions will at some level be limited by this 
shift. An intriguing alternative would be to instead use 
a fermionic atom, which should have no collisional shift, 
and thus density induced frequency noise will not be an 
issue. On the other hand, the lack of collisions may also 
lead to more rapid decoherence: collisions appear to pre- 
serve our bosonic system from the decohering effects of 
spatial inhomogeneity. 

We acknowledge useful conversations with the other 
members of the JILA BEC collaboration. This work was 
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